On the Enhanced Interstellar Scattering Toward B 1849+005 
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ABSTRACT 

This paper reports new Very Large Array (VLA) and Very Long Baseline Array (VLBA) 
observations of the extragalactic source B1849+005 (Galactic coordinates, £ = 33? 44, b = +0?21) 
at frequencies between 0.33 and 15 GHz and the re-analysis of archival VLA observations at 0.33, 
1.5, and 4.9 GHz. The structure of this source is complex, confirming previous suggestions, but 
interstellar scattering dominates the structure of the central component at least to 15 GHz. An 
analysis of the phase structure functions of the interferometric visibilities shows the density fluc- 
tuations along this line of sight to be anisotropic (axial ratio = 1.3) with a frequency-independent 
position angle, and having an inner scale of roughly a few hundred kilometers. The anisotropics 
occur on length scales of order 10 15 cm (D/5 kpc), which within the context of certain magneto- 
hydrodynamic turbulence theories indicates the length scale on which the kinetic and magnetic 
energy densities are comparable. A conservative upper limit on the velocity of the scattering 
material is 1800 km s _1 , based on the lack of changes in the shapes of the 0.33 GHz images. 
In the 0.33 GHz field of view, there are a number of other sources that might also be heavily 
scattered, which would suggest that there are large changes in the scattering strength on lines of 
sight separated by a degree or less. Both B1849+005 and PSR B1849+00 are highly scattered, 
and they are separated by only 13'. If the lines of sight are affected by the same "clump" of scat- 
tering material, it must be at least 2.3 kpc distant. However, a detailed attempt to account for 
the scattering observables toward these sources (angular broadening of the extragalactic source, 
pulse broadening of the pulsar, and upper limits on the angular broadening of the pulsar) does 
not produce a self-consistent set of parameters for a clump that can reproduce all three mea- 
sured scattering observables. A clump of Ha emission, possibly associated with the H II region 
G33. 418— 0.004, lies between these two lines of sight, but it seems unable to account for all of the 
required excess scattering. 



Subject headings: 
(PSR B1849+00) 

Introduction 



galaxies: individual (B1849+005) - - ISM: structure 
radio continuum: general — scattering 



pulsars: individual 



Density fluctuations on subparsec scales in the 
interstellar plasma produce refractive index fluc- 
tuations. In turn these refractive index fluctua- 
tions scatter radio radiation propagating through 
the plasma. The result is a rich set of observables 
including angular broadening of compact sources, 
pulsar intensity scintillations, and low- frequency 
variability of extragalactic sources (for reviews see 
Rickett 1990; Narayan 1992). 

The density fluctuations responsible for inter- 



stellar scattering appear to be well parameter- 
ized by a power-law power spectrum. In many 
cases the spectral index of the power spectrum 
is near 11/3, the value predicted by the Kol- 
mogorov theory of turbulence for neutral, incom- 
pressible unmagnetized fluids. As the interstel- 
lar plasma is compressible, ionized, and magne- 
tized, it would appear to violate these assump- 
tions. There has been recent progress in under- 
standing how a Kolmogorov-like spectrum could 
arise in a magnetized plasma, though (Lithwick & 
Goldrcich 2003, and references within). Combined 



with the large Reynolds number of the interstellar 
plasma, this suggests that the density fluctuations 
arise from a turbulent process. A small number of 
lines of sight toward pulsars are notable counter- 
examples, having spectral indices larger than 4 
(Hcwish, Wolszczan, & Graham 1985; Wolzszcan 
& Cordes 1987; Gupta, Rickett, & Lyne 1988; 
Clegg, Fiedler, & Cordes 1993) or an index ap- 
proaching 4 (Spangler & Gwinn 1990; Moran et al. 
1990; Wilkinson, Narayan, & Spencer 1994; Mol- 
nar et al. 1995). An index larger than 4 arises 
from a random distribution of density fluctuations, 
while an index approaching 4 is interpreted as a 
modification of the power-law spectrum by an in- 
ner scale indicative of a dissipative scale for the 
turbulence. Thus, determination of the spectral 
index of the density fluctuation spectrum can pro- 
vide clues about the mechanisms by which the 
density fluctuations originate. 

Wilkinson et al. (1994) showed that not only 
can a density spectral index be obtained, but that 
multi-frequency observations may probe the orga- 
nization and orientation of the density fluctuations 
responsible for scattering. They showed that the 
axial ratio and position angle of Cyg X-3 are fre- 
quency dependent, a frequency dependence that 
they interpreted in terms of an increasingly or- 
dered magnetic field on smaller and smaller scales. 
Trotter, Moran, & Rodriguez (1998) showed sub- 
sequently that the axial ratio and orientation of 
the scattering disk of NGC 6334B also exhibited 
a frequency dependence. Thus, scattering obser- 
vations may be a probe of the magnetic field on 
subparsec scales. 

Conversely, OH masers seen toward the Galac- 
tic center exhibit some of the most extreme 
anisotropic images seen anywhere in the Galaxy 
(van Langevelde et al. 1992; Frail et al. 1994). 
While it is not possible to measure a frequency 
dependence for the scattering disks of the OH 
masers, the scattering disk of Sgr A* exhibits no 
change in its axial ratio or orientation over nearly 
1.5 decades in frequency (Lo et al. 1985; Jauncey 
et al. 1989; Alberdi et al. 1993; Yusef-Zadeh et al. 
1994; Marcaide et al. 1999; Reid et al. 1999). It is 
not yet clear what the different frequency behav- 
ior for the scattering diameters of Cyg X-3 and 
NGC 6334B vs. Sgr A* represents. Possibilities 
might include a difference between the scattering 
media in the Galactic center and disk or a much 



stronger magnetic field in the Galactic center. 

The source B1849+005 (Galactic coordinates, 
£ = 33?44, b = +0?21) is seen along the third 
most highly scattered line of sight known, after 
the Galactic center (Lazio & Cordes 1998) and 
NGC 6334B (Trotter et al. 1998). Spangler et al. 
(1986) conclude that it is extragalactic because it 
is a compact, flat-spectrum variable source. Fur- 
ther evidence that B1849+005 is extragalactic is 
found by comparing its H I absorption spectrum 
to that of the pulsar PSR B1849+00, which is 13' 
from B1849+005. It must be more distant than 
PSR B1849+00 and is probably at least 25 kpc 
distant (van Gorkom et al. 1982; Dickey et al. 
1983; Clifton et al. 1988; Frail & Weisberg 1990). 
PSR B1849+00 is also the most heavily scat- 
tered pulsar known, with a pulse broadening time 
of 0.2 seconds at 1.4 GHz. Thus, comparison of 
the line of sight to B1849+005 and PSR B1849+00 
can probe changes in the scattering medium trans- 
verse to the line of sight on scales not normally 
accessible. 

This paper reports multi-frequency observa- 
tions of B1849+005, spanning nearly two decades 
in frequency, that are used to reassess the struc- 
ture and nature of B1849+005 and constrain the 
scattering properties along this line of sight. In 
§2 I describe the observations and data reduction, 
in §3 I use these observations to extract various 
scattering parameters along the line of sight to 
B1849+005, in §4 I combine my observations with 
those from the literature to probe the scattering 
environment around B 1849+005 and particularly 
toward PSR B1849+00, and in §5 I summarize my 
results. 

2. Observations and Analysis 

The scattering disk of B1849+005 is large 
enough that it can be resolved by a range of fre- 
quencies and baselines accessible with the VLA 
and VLBA. The details of the observations and 
data reduction for the VLA and VLBA observa- 
tions are sufficiently different that I discuss them 
separately. Table 1 presents the observing log. 



Table 1 
Observing Log 
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2.1. VLA Observations (0.33 4.9 GHz) 

2.1.1. 0.33 GHz 

Observations were obtained at 0.33 GHz with 
the A-configuration VLA on 1998 March 19. The 
sources 3C48 and 3C286 were used in calibrat- 
ing the visibility amplitudes, and the VLA cal- 
ibrator J1859+1259 was used in calibrating the 
visibility phases. Standard calibration and radio 
frequency interference (RFI) excision procedures 
were used within the NRAO Astronomical Image 
Processing System (aips). Several iterations of 
self-calibration and imaging were used to approach 
the thermal noise level. During imaging, a poly- 
hedron algorithm (Cornwell & Perley 1992) was 
used to account for the non-coplanar nature of the 
VLA. 

At 0.33 GHz the field of view is large (2? 5 
FWHM) and contains many sources. In order to 
isolate the visibilities from B1849+005, the follow- 
ing procedure was used. Sources within the pri- 
mary beam field of view as well as strong sources 
outside the primary beam were identified in the 
NRAO VLA Sky Survey (Condon ct al. 1998). 
Seventy three such sources were identified, in- 
cluding B1849+005. Following Perley (1989), the 
Fourier transform of all of the sources was sub- 
tracted from the visibility data, a modest amount 
of additional data editing was performed, and the 
Fourier transform of only B1849+005 was added 
back to the visibility data. 

The imaging was done with the visibility data 
at the spectral resolution at which it was acquired 
(97.7 kHz channels), both to avoid bandwidth 
smearing and to make RFI easier to identify. After 



removing all other sources, bandwidth smearing is 
no longer an issue, and the data were averaged 
over frequency. Also, in order to reduce the data 
volume before further processing, the visibilities 
were averaged from 10 s to 240 s. This time and 
frequency averaging also further "dilutes" the im- 
pact of any unsubtracted sources or errors in the 
models of the sources that were subtracted. 

I also obtained from the VLA Archive (pro- 
gram AS305) A-configuration VLA observations 
of B1849+005 acquired originally on 1987 Au- 
gust 20 by Fey, Spangler, & Cordes (1991). The 
calibration of these observations was similar (but 
with the VLA calibrator B1938— 155 used to cali- 
brate the phases), but, given the relatively limited 
amount of visibility data (15 min. observation), 
only B1849+005 was imaged, rather than using 
targeted facetting. The limited visibility data also 
implies a higher noise level and limited image fi- 
delity. 

Figure 1 shows the resulting 0.33 GHz images 
of B1849+005. 

2.1.2. 1.5 GHz 

The nominal resolution of the A-configuration 
VLA at 1.5 GHz is approximately l'/5. The 
expected scattering diameter of B1849+005 at 
this frequency is approximately 0'.'3. "Super- 
resolution" imaging can provide better than the 
nominal resolution by appropriate weighting of the 
visibility data at the cost of a modest reduction in 
sensitivity (Briggs et al. 1999). Given the strength 
of this source (~ 0.5 Jy at 1.5 GHz), this is an ac- 
ceptable trade-off. 

I obtained from the VLA Archive (program 



AS305) A-configuration VLA observations of 
B1849+005 acquired originally on 1987 August 20 
by Fey et al. (1991) with the aim of imaging the 
source at super-resolution. The source 3C 286 
was used to calibrate the visibility amplitude, and 
the VLA calibrator B1821+107 was used in cali- 
brating the visibility phases. Standard calibration 
procedures were used, including a number of iter- 
ations of phase self-calibration. 

As was done in imaging the 0.33 GHz observa- 
tions, a polyhedral algorithm was used to account 
for the non-coplanarity of the VLA and image the 
entire field of view (30' FWHM). In practice, the 
number of other sources in the field was sufficiently 
small and their flux densities sufficiently low that 
ignoring the VLA's non-coplanarity would have 
made little difference in the final image. 

Experiments with the amount of super- resolution 
used showed that reasonable choices for the 
amount of weighting led to resolution improve- 
ments by as much as 20% (i.e., decrease in beam 
diameter). Figure 2 shows the resulting image 
of B1849+005 with a modest amount of super- 
resolution (10% reduction in the beam diameter). 
In practice, the amount of superresolution used 
made little difference in the measured image diam- 
eter, but the measured image diameter is smaller 
than the beam diameter. Decreasing the beam 
diameter increases the ratio of the image diameter 
to the beam diameter and therefore increases the 
robustness of the measured image diameter, par- 
ticularly if the image diameter is smaller than the 
beam diameter, as is the case here. 

2.1.3. 4.9 GHz 

Fey et al. (1991) also obtained a single hour an- 
gle (1 min. duration) on B1849+005 at 4.9 GHz 
during the same set of observations used to acquire 
the 1.5 GHz observations. Based on the results of 
imaging the 1.5 GHz observations, these 4.9 GHz 
observations were re-analyzed also. The source 
3C 286 was used to calibrate the visibility ampli- 
tude, and the VLA calibrator B1821+107 was used 
in calibrating the visibility phases. Standard cal- 
ibration procedures were used, including a small 
number of iterations of phase self-calibration. Fig- 
ure 3 shows the resulting image. 
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Fig. 1.— B1849+005 at 0.33 GHz. Left The image 
resulting from the 1987 August 20 observations. 
The rms noise level is 80 mJy beam , and the 



contour levels are 80 mJy beam 



-3, 3, 5, 7.07, 



10, .... The beam is 15'.' 7 x 5'.'1 and is shown in 
the lower left. Right The image resulting from 
the 1998 March 19 observations. The rms noise 
level is 4.4 mJy beam -1 , and contour levels arc 
4.4 mJy beam" 1 x -3, 3, 5, 7.07, 10, .... The 
beam is 7'.' 4 x 4'.' 9 and is shown in the lower left. 
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Fig. 2.— B1849+005 at 1.5 GHz, imaged using 
a modest amount of super-resolution to obtain a 
10% reduction in the nominal beam diameter. The 
rms noise level is 0.32 mJy beam , and contour 
levels are 0.32 mJy beam -1 x -3, 5, 7.07, 10, 
14.1, 20, ... . The beam is l'/3 x l'/3 and is shown 
in the lower left. 



2.2. VLBA Observations (2.3-22 GHz) 

Observations at 5, 8.4, 15, and 22 GHz were 
obtained with the VLBA and one antenna of the 
VLA on 1998 March 29. Observations at 2.3 GHz 
were obtained with only the inner antennas of 
the VLBA on 2002 August 23. No fringes were 
detected at 22 GHz. The lack of a detection 
at 22 GHz is consistent with the lower sensitivity 
of the VLBA at this frequency relative to the other 
frequencies observed. At the other frequencies, 
the correlated visibilities were fringe-fit and am- 
plitude calibrated in the standard fashion within 
AlPS. Fringe-fitting intervals were 2.5-3 min. The 
calibrated visibilities were then exported to the 
Caltech dif map program for hybrid mapping. 

For the 5, 8.4, and 15 GHz observations, af- 
ter obtaining an image with a noise near, but 
still well above, the thermal limit, I subtracted 
the Fourier transform of B1849+005 from the vis- 
ibility data, performed a modest amount of ad- 
ditional data editing, and then added the Fourier 
transform of B1849+005 back to the visibility data 
(Perley 1989). Figures 4-6 show the resulting im- 
ages of B1849+005 between 5 and 15 GHz. 

At 2.3 GHz, a significant detection of B1849+005 
was obtained only on the shortest baseline, be- 
tween the LA and PT antennas. As shown be- 
low, this source has a fairly simple structure on 
milliarcsecond scales. In order to estimate a di- 
ameter, the correlated visibility amplitude on the 
LA-PT baseline was compared to total flux den- 
sity of the source. Fiirst et al. (1990) measure 
the flux density of the source at 2.7 GHz to be 
0.77 Jy; flux densities obtained by interpolating 
its spectrum between 1.5 and 5 GHz (Fey et al. 
1991; this work) range from 0.6 to 0.8 Jy. Because 
the measurement by Fiirst et al. (1990) utilized a 
single-dish telescope, I adopt 0.77 Jy as its flux 
density at 2.3 GHz. Assuming the source to be 
a circular Gaussian, its diameter is then 110 mas 
(Table 2). 

For further analysis, the data were averaged in 
time. The averaging times were 30 s at 2.3 GHz, 
60 s at 5 GHz, 30 s at 8.4 GHz, and 15 s at 15 GHz. 
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Fig. 3.— B1849+005 at 4.9 GHz, as observed by 
the VLA. The rms noise level is 1 mJy beam , 
and contour levels are 1 mJy beam -1 x —3, 5, 
7.07, 10, 14.1, .... The beam is 0'.'5 x 0'.'4 and is 
shown in the lower left. 
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Fig. 4. — The central component of B1849+005 
at 5 GHz, as observed by the VLBA. The rms 
noise level is 0.39 mJy beam -1 , and contour levels 
are 0.39 mJy beam -1 
The beam is 14 mas 
the lower left. 
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Fig. 5. — The central component of 

B1849+005 at 8.4 GHz. The rms noise level 
is 0.51 mjy beam -1 , and contour levels are 
0.51 mjy beam- 1 x -3, 3, 5, 7.07, 10, .... The 
beam is 4.2 mas x 2.8 mas and is shown in the 
lower left. 
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Fig. 6. — The central component of 

B1849+005 at 15 GHz. The rms noise level 
is 0.84 mjy beam -1 , and contour levels are 
0.84 mjy beam -1 x -2, 3, 5, 7.07, 10, .... The 
beam is 2.0 mas x 0.88 mas and is shown in the 
lower left. 



3. Interstellar Scattering Toward B 1849+005 
3.1. Relevant Scattering Formulae 

The intcrfcrometric visibility of an infinitely 
distant point source seen through a medium filled 
with random density fluctuations is 



V(b) 



-D*(b)/2 



(1) 



for an interferometer with baseline length b. For 
the density power spectrum relevant in the inter- 
stellar medium (Armstrong, Rickett, & Spanglcr 
1995), the form of the phase structure function 
depends upon the value of the importance of l\ 
relative to the typical baseline length b. Antici- 
pating later results, I consider two possibilities for 
D^b) (e.g., Coles et al. 1987; Cordes & Lazio 1991; 
Lambert & Rickett 1999) 

D^b) = 87rr;? A 2 SM/ (a) x 

b a - 2 b > I 

r(a/2) ir 2 [iF 1 (a/2; 1; -(b/h) 2 ) - 1] b < h 

Here r e is the classical electron radius, A is the 
free-space observing wavelength, /(a) is a slowly- 
varying function of a with a value near unity, 
\Fi(a,b;x) is the confluent hypergeometric func- 
tion, and the scattering measure SM is propor- 
tional to the linc-of-sight integrated rms electron 
density. (Minor modifications are required if the 
source is not infinitely distant, as assumed here.) 
For reference, in the local interstellar medium, 
SM - 10 -35 kpc m" 20 / 3 for a 1 kpc path length. 
The resulting scattering diameter, for a dis- 
tant source and assuming a = 11/3 (Rickett 1990; 
Cordes & Lazio 1991), is 



\2) 



128 mas SM 



0.6 ,,-2.2 
y GHz 



71 mas SM 1/2 ^ 2 Z ( 



h 



-1/6 



100 km; 



&> h, 
6< h. 

(3) 
where z^ghz is the observing frequency in GHz and 
SM is in its canonical units of kpc m~ 20 / 3 . 

3.2. Angular Broadening and the Appar- 
ent Structure of B1849+005 

On the basis of its spectrum, Fey et al. (1991) 
suggested that the structure of B 1849+005 may be 
more complex than that of a single compact com- 
ponent. Figures 2 and 3 demonstrate that this is 
the case. Although a strong central component 



does dominate the structure of the source, there 
are weaker components both to the southwest and 
northeast. The components to the northeast are 
at best only faintly visible in the VLA observa- 
tions at 4.9 GHz (Figure 3), though their absence 
probably results from the limited hour-angle cov- 
erage and surface brightness sensitivity of the A- 
configuration VLA at this frequency. The central 
component may show a hint of a jet emerging to 
the southeast and bending toward the southwest. 
The existence of a jet must be considered tenta- 
tive, though, given the extremely limited hour an- 
gle coverage (only 1 min. observation). As the 
0.33 GHz image is elongated in approximately the 
same orientation as that of the central and south- 
west components, I have also re-imaged the (1998 
March 19) 0.33 GHz observations utilizing super- 
resolution. There is no indication of this second 
component, though the orientation of the beam is 
similar to that of the central and southwest com- 
ponents. 

The expressions of §3.1 assume an infinitely dis- 
tant point source. To what extent can they be 
applied to B 1849+005, given its complex struc- 
ture? Figure 7 shows the diameter of the cen- 
tral component of B 1849+005 as a function of 
observing frequency. Although Figures 1-6 show 
images of B1849+005, diameters measured from 
these images were not used in constructing Fig- 
ure 7. Rather, Gaussian models were fit to the 
visibility data, except in the case of the 2.3 GHz 
observations for which a circular gaussian was as- 
sumed (§2). Also plotted in Figure 7 are mea- 
surements at various other frequencies. At these 
other frequencies (0.25 and 0.408 GHz) the visibil- 
ity data were not available for analysis (as below), 
but angular diameters could be determined. For 
completeness, Table 2 summarizes the various fits 
for the diameter of the source. 

Figure 7 also shows three curves, ex ^~ 2,2 , 
6 ex v~ 2 , and 8 oc v~ 1A . The first two are relevant 
for strongly scattered sources, with the first be- 
ing appropriate if the electron density spectrum 
is a power law with a Kolmogorov spectral in- 
dex, a = 11/3, while the second is appropriate 
if the range of observing baselines is comparable 
to l\. The third curve shows the frequency de- 
pendence for the angular broadening if scattering 
toward B1849+005 is anomalous in a manner sim- 
ilar to that seen for the pulsar PSR B1849+00 
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Fig. 7. — The angular diameter of B1849+005 as 
a function of observing frequency. The angular di- 
ameter was obtained by fitting a Gaussian either 
to the image or to the visibility; if an elliptical 
Gaussian was used in the fit, shown is the geomet- 
ric mean of the major and minor axes. At most fre- 
quencies, the formal uncertainties are comparable 
to the size of the points. The datum at 0.25 GHz 
is from P. Rao (2001, private communication) and 
that at 0.408 from Dennison ct al. (1984). For ref- 
erence, three curves are also plotted, one showing 
a v~ 2 dependence for the diameter (dotted line), 
one showing a v~ 22 dependence (dashed line), and 
one showing a v~ 1A dependence (dot-dash line) 
that would be expected if the angular broadening 
for B 1849+005 is anomalous in a manner similar 
to that for PSR B1849+00 (§4.2). All three curves 
are constrained to go through the diameter mea- 
sured at 5 GHz. 



(§4.2). The relatively good agreement between 
the angular diameters and the first two curves is 
an indication that scattering dominates the ap- 
parent structure of the central component at least 
up to 15 GHz. Interpolating the angular diame- 
ter to the fiducial frequency of 1 GHz, the scat- 
tering diameter of B1849+005 is 0'.'6. Based on 
their more limited set of observations, Fey et al. 
(1991) reached a similar conclusion: B1849+005 
is heavily scattered, though exhibits a more com- 
plex structure than simply a single compact com- 
ponent. 

In addition, the angular diameter of the south- 
west component was measured, using the images 
shown in Figures 2 and 3. Its diameter is approx- 
imately constant between 1.5 and 4.9 GHz, indi- 



Table 2 
Gaussian Model Fits to the Central Component of B1849+005 



Frequency 


S 


a 


b 


$ 


Reference 


(GHz) 


(Jy) 


(mas) 


(mas) 


(°) 




0.25 


1.624 ± 0.009 


11 400 ± 70 


8430 ± 60 


44.2 ± 0.7 


2 


0.33 


2.1 ± 0.2 


6300 ± 1000 


2080 ± 100 


18 ± 7 


1 


0.33 


2.12 ± 0.02 


7050 ± 350 


2560 ± 100 


25.8 ± 0.7 


1 


0.41 


1.7 ± 0.2 


> 600 






3 


1.5 a 


0.63 ± 0.06 


440 ± 6 


190 ± 6 


161 ± 1 


1 


2.3 


0.77 ± 0.1 


110 ± 5 


110 ± 5 




1 


5.0 


0.739 ± 0.009 


19.34 ± 1 


17.32 ± 1 


27.0 ± 0.6 


1 


8.1 


0.634 ± 0.006 


6.71 ± 0.3 


4.91 ± 0.2 


-33. 2± 0.3 


1 


15 


0.831 ± 0.008 


2.017 ± 0.1 


1.623 ± 0.08 


-21. 1± 0.5 


1 



a Only the central component was fit. 

Note. — The visibility data are fit with a single Gaussian model having a flux density S, 
major axis a, minor axis b, and position angle <1>. In some cases a circular Gaussian was 
fit; no position angle is reported for these cases. The first entry for 0.33 GHz is for the 
1987 August 20 observations; the second is for the 1998 March 19 observations. The stated 
uncertainties are the formal statistical uncertainties for the fits; the actual uncertainties due 
to systematic effects are almost certainly larger. 

References. — (1) This work; (2) A. P. Rao (2002, private communication); (3) Dennison 
ct al. (1984) 



eating that its diameter is much larger than the 
scattering diameter. This is consistent with its 
somewhat diffuse appearance at 4.9 GHz. 

3.3. The Electron Density Spectrum 

A more profitable analysis of the observations 
is to use the structure function (equations 1-2) to 
extract parameters of the electron density spec- 
trum. My analysis follows closely that of Spanglcr 
& Cordes (1998). Figure 8 shows an example of a 
structure function constructed from the 8.4 GHz 
visibility data. 

Clearly apparent in Figure 8 is a non-random 
scatter in the structure function, e.g., around 
10 MA. Moreover, Figures 1 and 4-6 show 
that the scattered image has (varying) degrees 
of anisotropy. This anistropy arises presumably 
because the density fluctuations responsible for 
the scattering are themselves anisotropic. 

The anisotropy of the image and visibility data 
(and of the underlying density fluctuations) are 
described by an orientation ^ and aspect ratio r\. 
In fitting the visibility data to equation (1), the 
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Fig. 8. — Structure function for the 8.4 GHz vis- 
ibility data using the isotropic baseline. A small 
number of visibilities below 2 MA are not shown. 

baseline b is replaced by the rotundate baseline 

b 2 r = (ucos* + usin*) 2 + ry 2 (-usin* + wcos*) 2 

(4) 
where u and v are the traditional interferome- 
ter baseline coordinates measured in units of the 
observing wavelength. The rotundate baseline is 



simply the baseline measured in an anisotropic co- 
ordinate system that is rotated relative to that of 
the interferometer's coordinate system. In this ro- 
tundate coordinate system, ^ is measured coun- 
terclockwise from the u axis, in contrast to the 
position image in an image, which is measured 
counterclockwise from north. Thus, there is a 90° 
offset from an image position angle and \l/. 

The non-random scatter in Figure 8 arises be- 
cause of the use of the isotropic baseline (rj = 1 
and ^f = 0°). Figure 9 illustrates the use of the ro- 
tundate baseline in plotting the structure function 
(using the best-fitting power-law density spectrum 
model found below). Clearly apparent is the dra- 
matic reduction in the scatter of the data; similar 
improvements are seen at the other frequencies. 
Figure 9 also demonstrates that the structure of 
B1849+005 does not have an adverse impact on ef- 
forts to extract information about the scattering. 
A single source component does an excellent job of 
fitting the visibilities. As might be expected from 
Figures 2 and 3, the VLBI observations are insen- 
sitive to the relatively low brightness temperature 
components to the southwest and northeast. 




10 20 

Rotundote Baseline (MX) 

Fig. 9. — Structure function for the 8.4 GHz visi- 
bility data using the rotundate baseline. The solid 
line shows the best fitting model for a power-law 
density spectrum (viz. Table 3). A small number 
of visibilities below 2 MA are not shown. 

One important aspect of the use of the ro- 
tundate baseline is that the structure functions 
in equation (2) have been derived assuming an 
isotropic density spectrum. This ex post facto re- 
placement of the rotundate baseline for the base- 
line violates that assumption. This replacement 



should be approximately valid provided that 77 is 
not much larger than unity. Moreover, Chandran 
& Backer (2002) have considered the propagation 
of radio waves through an anisotropic medium. 
They find that expressions like those of equa- 
tion (2) continue to hold; in particular, the power- 
law dependences remain unchanged. 

Using rotundate baselines and the structure 
functions of equation (2), there are four or five pa- 
rameters for which to fit, depending upon the rel- 
ative importance of the inner scale — a, SM, 77, 'J, 
and l\. I used a genetic algorithm (Lazio 1997) to 
find the best fitting parameters in a (unweighted) 
minimum chi-squared sense at each frequency. 

I have fit the structure functions (as in Figure 8) 
for the VLBI observations (2.3-15 GHz) and the 
0.33 GHz VLA observations of 1998 March 19. 
The visibility data from the 1.5 or 4.9 GHz ob- 
servations from the VLA were not used because it 
is clear that the source structure is complex and, 
at 4.9 GHz, the expected scattering diameter is 
much smaller than what can be resolved by the 
VLA. Also, the 0.33 GHz visibility data from the 
1987 August 20 VLA observations were not used 
because of their limited quantity and restricted u- 
v coverage. 

Not all of the data were used in the fitting. 
At long baselines, the structure functions saturate 
while at short baselines they reach noise plateaus 
(viz. equations 1 and 2). The saturation at long 
baselines is visible near 30 MA in Figures 8 and 9; 
the noise plateau at short baselines occurs be- 
low 2 MA (for the 8.4 GHz data). Consequently, 
limits were set at both short and long baselines. 
After excluding noisy data, the total number of 
data available to fit ranged from several hundred 
to several thousand. 

Figures 10 and 11 show the best fitting param- 
eters and their confidence ranges as a function of 
frequency for both a power-law density spectrum 
(b ^ l\) and a density spectrum with an inner 
scale (b < l\), respectively; Tables 3 and 4 tabu- 
late the best-fitting model parameters. Figures 10 
and 11 also indicate the (weighted) mean values, 
over all frequencies, for the various parameters. A 
X 2 statistic was used to test whether the best fit- 
values are consistent with a single mean value, i.e., 
no frequency dependence. 
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Fig. 10. — The best fit parameters as a function of frequency for a power-law density spectrum. The 
confidence regions shown are the amount by which the parameter must change in order to cause the chi- 
squared to increase by unity from its minimum value. The horizontal dotted line shows the weighted 
mean value for the parameter. Because B1849+005 was detected on only a single baseline at 2.3 GHz, the 
anisotropy i] and its orientation ^> were constrained to be 1 and (cos 'I' = 1), respectively. At 1.5 GHz, a 
structure function analysis was not conducted; the values shown are derived from analysis of the image. 
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Fig. 11. — The best fit parameters as a function of frequency for a density spectrum with an inner scale. As 
for Figure 10. 
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The confidence ranges for the parameters have 
been estimated both by a local grid search and a 
non-linear least squares minimization around the 
minimum \ 2 ■ I n the former method, the confi- 
dence range for a parameter was determined by 
varying that parameter, while holding all others 
fixed at their best fitting values, until the \ 2 ap- 
proximately doubled from its minimum value. In 
the latter method a Levenberg-Marquardt mini- 
mization routine was used to determine the covari- 
ance matrix of the parameters. The Levenberg- 
Marquardt minimization yields confidence ranges 
much smaller than the grid search. We have 
adopted the confidence ranges for the grid search; 
adopting smaller confidence ranges would make 
our following conclusions more robust. 

One aspect of Figures 10 and 11 warrants com- 
ment. The various parameters — a, SM, 77, "f, 
and l\ — are plotted as functions of frequency 
These parameters were introduced as frequency- 
independent quantities (equation 2). For instance, 
SM is proportional to the rms electron density. 
The apparent frequency dependence arises from 
the use of interferometers whose baseline lengths 
were not scaled with frequency. For an interferom- 
eter with a typical (and fixed) baseline length b, it 
will resolve the scattering disk only if the scatter- 
ing disk is comparable in size to the typical angle 
measured by the interferometer X/b. In turn if the 
scattering medium is at a distance D, the scatter- 
ing disk will be produced by density fluctuations 
on scales of/' ~ D(X/b). This length scale is man- 
ifestly frequency dependent. Thus, an apparent 
frequency dependence in a, SM, 77, ^>, or li would 
reflect different parameter values being obtained 
on different length scales within the spectrum of 
scattering irregularities. For brevity, in the con- 
text of Figures 10 and 11,1 shall continue to speak 
of a frequency dependence for a, SM, 77, 'J, and l\ 
with the understanding that any frequency depen- 
dence arises from the different length scales probed 
in the scattering medium. 

3.3.1. Image Anisotropy 

Regardless of whether a power-law or inner- 
scale density spectrum is fit to the structure func- 
tions, the result confirms what the casual inspec- 
tion of Figures 1-6 suggest: The scattering toward 
B1849+005 is decidedly anisotropic. 

At 0.33 GHz, the anisotropy appears large, 



77 ~ 2.3, but the weaker components, especially 
the one to the southest, may be contaminating the 
analysis. The beam diameter at 0.33 GHz is com- 
parable to the separation between the components 
and the orientation of the source at 1.5 GHz is 
similar to that of both the source and the position 
angle of the beam at 0.33 GHz. The visibilities 
at 0.33 GHz appear consistent with only a single 
component being present (§3.2 and Figure 7), but 
the presence of the other components cannot be 
excluded. Indeed, convolving the 1.5 GHz image, 
Figure 2, to a resolution comparable to that of 
the 0.33 GHz images, Figure 1, produces an im- 
age similar to the 0.33 GHz images. 

Even if the analysis at 0.33 GHz is disregarded 
(and at 1.5 and 2.3 GHz, for which less compre- 
hensive analyses could be performed), the higher 
frequency data show a clear anisotropy. The fits 
at 5.0 and 15 GHz are just marginally consistent 
with an isotropic image, while the fit at 8.4 GHz 
excludes an isotropic image. The mean value for 
the image elongation is fj = 1.2. This value is 
comparable to that measured for Cyg X-3 (« 1.3, 
Wilkinson et al. 1994; Molnar et al. 1995) and 
NGC 6334B (1.2-1.5, Trotter et al. 1998), within 
the range of anisotropics measured for the extra- 
galactic sources seen through the Cygnus region 
(1.1-1.8, Spangler & Cordes 1988, 1998; Desai & 
Fey 2001), and somewhat smaller than the value 
measured for the OH masers toward W49N (2-3, 
Desai, Gwinn, & Diamond 1994) and the Galactic 
center (« 2.5, Frail et al. 1994). Romani, Narayan, 
& Blandford (1986) have predicted the anisotropy 
that would be induced, even if the small-scale 
(diffractive) density fluctuations responsible for 
angular broadening are isotropic, by refractive ef- 
fects from large-scale density fluctuations. For a 
thin scattering screen with a Kolmogorov density 
spectrum, the expected anisotropy is (77) < 1.05. 
The scattering toward B1849+005 is both strong 
and anisotropic, and the anisotropy must result 
from the density fluctuations responsible for the 
scattering. 

Whether the position angle ^ of the anistropy 
changes with frequency, as Wilkinson et al. (1994) 
and Trotter et al. (1998) found for Cyg X-3 and 
NGC 6334B, respectively, is less certain. In both 
Figure 10 and 11, an apparent trend of a decreas- 
ing position angle (or increasing cos , f) is seen. 
This is the case regardless of whether the 0.33 GHz 
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Table 3 
Power-law Density Spectrum Model Fitting 



Frequency 
(GHz) 


Q 


V 


(°) 


logSM 
(kpc m- 20 / 3 ) 


0.33 

1.5 a 


3.68JS3? 


2.29+°f 
2.3 


1581^4 
161 


47+ - 41 

"•^'-0.06 

0.98 


2.3 b 

5.0 
8.4 
15 


3.60 

o 04+O.OI 
oq+0.02 

°-°°-o.io 

01+001 

°- yi -0.03 


1 10+0.O6 
J-^-o.is 

1 4Q+0-08 

1 29+ 007 


134+46 
119+J? 

114+2 y 


0-60+8:8 

0-79+°;°« 

r, 7 q + 0.07 

u -'°-0.21 


mean 


3.92 


1.30 


141 


0.68 



a At this frequency, the scattering parameters were derived 
from the image, not from a structure function analysis. 

b At this frequency, the scattering diameter was assumed to 
be circular. 



Table 4 
Power-law with Inner Scale Density Spectrum Model Fitting 



Frequency 


a 


'/ 


* 


logSM 


h 


(GHz) 






(°) 


(kpc m- 20 / 3 ) 


(km) 


0.33 


7n +o.oi 

°-' u -0.10 


2.00+™ 5 


160+_ 4 7 


0.92l°;g 


190±£° 


2.3 


3.60 






1.30+0;! 


62_ 12 


5.0 
8.4 
15 


3.80±g;S 

q 7n +0.08 
°-' u -0.06 

3.80+°;° 


1 20+ 010 

i - zu -0.20 

1 sn+°- 30 

i.ou_ 20 
1 20+ - 04 


142+1 

125±g 

-191+32 
izi_ 18 


1 14+ - 21 
1 03+ 025 

±.UO_ 25 

1 16+ - 04 

J — LO -0.58 


500t° 00 
470+ 3 ° 8 
490+™ 


mean 


3.78 


1.23 


148 


1.14 


470 



*At this frequency, the scattering diameter was assumed to be circular. 
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datum is included or not. However, as both figures 
show, the uncertainties are also consistent with no 
change in position angle. 

Using only the higher frequency data, the scat- 
tering diameter of B1849+005 ranges from 19 mas 
at 5.0 GHz to 2 mas at 15 GHz. If the ma- 
terial responsible for the scattering is at a dis- 
tance of 5 kpc (§4.2), the size of the region on 
which the anisotropic density fluctuations occur 
must range from 1.5 x 10 14 -1.4 x 10 15 cm (10- 
95 AU) . These scales are comparable to those seen 
for NGC 6334B, though this is due partially to a 
selection effect in that the same telescope was used 
in both this work and that of Trotter et al. (1998) 
so the range of accessible spatial scales is the same. 

Within the context of the models for magne- 
tohydrodynamic turbulence developed by Goldrc- 
ich and collaborators (e.g., Lithwick & Goldrc- 
ich 2001), these length scales would represent the 
lower limit to their "MHD scale." On this length 
scale, the kinetic energy density (ex v 2 ) of the tur- 
bulent eddies has diminished to the point that the 
magnetic energy density (oc B 2 ) is comparable to 
it. Thus, this length scale marks where the tur- 
bulence makes a transition from isotropic on large 
scale (comparable to the outer scale) to increas- 
ingly anisotropic on small scales. 

3.3.2. Inner Scale 

Comparison of Figure 10 and 11 shows that the 
density spectral index a appears to change be- 
tween 0.33 GHz and 5.0 GHz for the power-law 
spectrum. At 0.33 GHz, the fit for a power-law 
spectrum finds a — 3.68 while above 5.0 GHz, 
a > 3.8. In contrast, the fit for an inner scale 
spectrum reveals a much more consistent set of 
values for a, with 3.7 < a < 3.8. Moreover, the 
fits for the inner scale itself find consistent values 
of a few hundred kilometers (though with large 
uncertainties). 

I regard this as a detection of the effects of the 
inner scale on scales of a few hundred kilometers. 
The maximum baselines of the VLA, used for the 
0.33 GHz observation, are 35 km. In contrast, the 
VLBA baselines, used for the higher frequency ob- 
servations, span the range from a few hundred to 
a few thousand kilometers. Thus the "break" be- 
tween 0.33 GHz and 5.0 GHz in the fits for a in the 
power-law density spectrum (Figure 10) reflects 



the range of baselines used in the observations. 

This value for the inner scale — a few hundred 
kilometers — is comparable to that found toward 
other sources (Spanglcr & Gwinn 1990; Moran et 
al. 1990; Wilkinson et al. 1994; Molnar et al. 1995). 

3.4. Secular Change of the Scattered Im- 
age 

If the plasma responsible for the scattering has 
a bulk (transverse) velocity v (and it remains 
relatively unchanged as it moves across the line 
of sight, "frozen- flow" approximation), changes 
in the image appearance should appear on time 
scales comparable to the refractive time scale t r = 
l r /v. The two observations at 0.33 GHz (Figure 1) 
provide a 10.6 yr baseline over which to search for 
secular changes in the image. 

Table 2 shows that both the fitted diameters 
and position angles for the source do not change 
appreciably. From both sets of observations 0d — 
4" at 0.33 GHz. The refractive scale is therefore 
6 x 10 16 cmDk P c, where £>kpc is the distance to 
the screen in kiloparsecs. There is little change 
in the image over the 10.6 yr interval. The ve- 
locity of the scattering material is therefore v < 
1800 km s~ -Dkpc- (Strictly, this velocity is a com- 
bination of the velocities of the source, medium, 
and observer. Here the velocity of the source is 
assumed to be effectively zero.) 

This upper limit is conservative. If the actual 
velocity were within a factor of two, say, of this 
upper limit, some modest changes in at least the 
orientation should start to become apparent. 

The source B1849+005 is only 14' from the su- 
pernova remnant (SNR) G33.6+0.1. Spangler et 
al. (1986) observed B1849+005 with the aim of as- 
sessing to what extent the SNR could contribute 
to the scattering of B 1849+005. While they in- 
ferred correctly that B 1849+005 is heavily scat- 
tered, they were reluctant to attribute all of the 
scattering to G33. 6+0.1 given the low latitude line 
of sight to B 1849+005. Unfortunately, this up- 
per limit on the velocity of the scattering material 
sheds no light on whether material from the SNR 
is responsible for the scattering. The large scat- 
tering diameter at 0.33 GHz means that relatively 
large length scales are probed within the scatter- 
ing material and that correspondingly long time 
intervals are needed to place stringent constraints 
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on the velocity. 

The higher frequency VLBI observations may 
prove useful in constraining the velocity of the 
scattering material in the future. For instance, 
repeating the 8.4 GHz observations after a 7-yr 
interval would place a 5 km s _1 limit on the ve- 
locity (presuming that no change in the scattering 
disk was seen). 

4. The Scattering Environment Around 
B1849+005 

4.1. Nearby Lines of Sight 

By virtue of the VLA's large field of view 
at 0.33 GHz (w 2?5), other highly scattered 
sources may be detectable in the direction of 
B1849+005. With only a single frequency, I 
cannot demonstrate unambiguously that scatter- 
ing dominates the structure of other sources in 
the field. However, given that B1849+005 and 
PSR B1849+00 (see following section) are both 
highly scattered and that PSR B1845-01 also 
exhibits a large amount of pulse broadening (Ra- 
machandran et al. 1997), it is reasonable to con- 
sider that scattering may contribute to the ob- 
served diameters of many of these sources. 

Figure 12 shows the angular diameters for all 
sources in the 0.33 GHz field of view of B1849+005 
that are judged to be compact and likely domi- 
nated by scattering. Table 5 tabulates the diame- 
ters of these sources. 
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Fig. 12.— Lines of sight in the 0.33 GHz field 
of view of B1849+005 that may also be dom- 
inated by scattering. The gray scale shows 
the smoothed, continuum-subtracted Ha emission 
from the Southern Ha Sky Survey Atlas (Gaustad 
et al. 2001) on a linear scale from to 30 Ralcighs. 
Sources marked with crosses appear unresolved; 
sources marked by circles appear resolved, with 
simple structure indicating that scattering proba- 
bly dominates; and sources marked by triangles in- 
dicate sources whose structures may be more com- 
plex and therefore contain contributions from both 
scattering and intrinsic structure. The size of a cir- 
cle or triangle is proportional to the source's diam- 
eter at 0.33 GHz. Sources identified are the pulsars 
PSR B1849+00; PSR B1854+00; PSR B1855-01; 
the supernova remnants Kes 79, 3C 391, and 
3C 392; and the mascr OH 34.26+0.15. The source 
B1849+005 is just above PSR B1849+00 (viz. Fig- 
ure 13). 
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Table 5 
Other Possibly Scattered Sources Near B1849+005 







Galactic 


Galactic 




Right Ascension 


Declination 


Longitude 


Latitude 


Diameter Structure 


(J2000) 






(") 


18 35 10.617 


01 45 52.86 


32.650 


4.422 


? 


18 37 18.258 


-01 50 26.68 


29.680 


2.302 


3.8 


18 41 44.789 


-01 52 52.35 


30.152 


1.297 


5.1 


18 43 58.583 


01 53 27.68 


33.766 


2.522 


3.7 


18 44 07.675 


00 33 17.10 


32.592 


1.879 


0.6 


18 44 50.864 


01 37 15.89 


33.624 


2.205 


3.4 


18 45 17.789 


00 07 42.94 


32.117 


1.308 


5.2 ? 


18 46 03.708 


00 03 39.24 


32.265 


1.168 


3.3 


18 46 05.681 


00 25 38.53 


32.703 


1.384 


3.1 


18 46 14.538 


-01 19 40.83 


31.158 


0.550 


5.0 


18 46 23.171 


00 18 28.74 


32.630 


1.264 




18 46 31.286 


00 09 08.35 


32.507 


1.163 




18 46 39.510 


00 16 23.36 


32.630 


1.188 




18 46 41.266 


00 55 14.98 


31.571 


0.637 




18 46 53.543 


00 41 02.80 


33.023 


1.323 


2.6 ? 


18 47 09.999 


01 41 55.15 


33.957 


1.724 




18 47 43.843 


01 54 33.63 


34.209 


1.695 


3.0 


18 47 55.076 


01 22 19.10 


33.752 


1.409 


7.0 


18 48 23.488 


-01 23 58.55 


31.339 


0.039 


4.2 


18 48 25.281 


00 46 36.09 


31.897 


0.317 


5.7 


18 48 34.187 


00 36 18.25 


32.067 


0.362 


5.5 


18 48 45.681 


01 52 55.68 


34.302 


1.453 


3.1 


18 49 12.392 


02 11 32.05 


34.629 


1.496 


4.1 


18 49 33.461 


01 54 29.11 


34.416 


1.288 


2.6 


18 50 10.398 


00 20 04.02 


33.086 


0.434 


2.7 


18 50 20.140 


00 49 13.91 


33.537 


0.619 


2.6 


18 50 40.088 


03 57 33.61 


36.370 


1.975 


2.4 


18 51 52.819 


00 40 01.42 


32.389 


-0.403 


3.9 


18 52 12.978 


02 22 52.19 


35.140 


0.912 


5.1 ? 


18 52 27.364 


00 32 00.70 


33.523 


0.017 




18 53 23.434 


00 21 37.42 


32.834 


-0.599 


1.5 


18 53 36.039 


00 27 14.06 


32.775 


-0.688 


1.9 


18 53 37.212 


02 24 56.68 


35.331 


0.616 


6.0 ? 


18 53 43.689 


00 19 56.05 


33.489 


-0.358 




18 54 34.997 


02 27 51.16 


35.484 


0.424 
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If the structures of these sources are dominated 
by scattering, there is substantial variation in the 
scattering on sub-degree angular scales. Slysh et 
al. (2001) have reached a similar conclusion by 
comparing their space VLBI observations of the 
maser OH 34.26+0.15 to scattering measurements 
for PSR B1849+00 and B1849+005. They find 
the maser to have an angular diameter of less than 
1 mas, from which they deduce a scattering mea- 
sure of 1.68 x 10~ 3 kpc m~ 20 / 3 , far less than that 
of PSR B1849+00 or B1849+005, though in this 
case it is likely that the maser lies in front of the 
strong scattering region. 

Figure 12 suggests that the region of intense, 
yet spatially-variable scattering may subtend a 
large angle on the sky. For instance, there are a 
number of sources near (18 h 47 I . n 5, +2°) whose 
diameters range from unresolved to 7". Nu- 
merous other cases of unresolved sources fairly 
close to sources with measurable diameters also 
can be identified. These lines of sight may be 
similar to other directions showing both heavy 
and anisotropic scattering and large changes in 
the scattering strength on degree-size scales such 
as the Cygnus region (Fey, Spangler, & Mutel 
1989; Lazio, Spangler, & Cordes 1990; Spangler 
& Cordes 1998; Desai & Fey 2001), the Galactic 
center (van Langevelde et al. 1992; Frail et al. 
1994), and along the Galactic plane Dennison et 
al. (1984); Fey et al. (1991). 

4.2. The Line of Sight to PSR B1849+00 

The pulsar PSR B1849+00 is only 13' from 
B1849+005 and is the most heavily scattered pul- 
sar known, with a pulse broadening time of 0.2 sec- 
onds at 1.4 GHz. Lohmer et al. (2001) show that 
between 1.4 and 2.7 GHz the pulse broadening of 
this pulsar is anomalous in the sense that its fre- 
quency dependence is v~ x with x = 2. 81^6 rather 
than the expected x « 4. In turn, this anomalous 
frequency scaling implies a density spectral index 
of a = 6.6^2'3- (For a review of other observations 
of this pulsar, see Stanimirovic et al. 2003.) 

My observations of B 1849+005 are at frequen- 
cies that bracket those used by Lohmer et al. 
(2001) to observe PSR B1849+00 (Figure 7). Un- 
fortunately the large uncertainties on the angu- 
lar diameters between 0.33 and 2.3 GHz pre- 
clude definite conclusions about the possibility 
of similar anomalous scattering along the line 



of sight to B1849+005. On the one hand, the 
frequency scaling for the angular diameter be- 
tween 1.5 and 2.3 GHz is v~ 2Z and the frequency 
scaling between 0.33 GHz and 5 GHz is v~ 2 . Both 
of these frequency scalings (as well as Figure 7) 
would suggest a frequency dependence of approx- 
imately v~ 2 is a much better description for the 
angular broadening than a dependence of approx- 
imately v~ 1A (because the pulse broadening and 
angular broadening are related as r ~ 8 2 ). Recall, 
though, that the 0.33 GHz diameter may be over- 
estimated due to contamination from the source 
structure (§3.2), the 1.5 GHz diameter is some- 
what smaller than the resolution of the array, and 
at 2.3 GHz the diameter has been estimated from 
only a single baseline. 

Cordes & Lazio (2001, their Figures 2 and 3) 
showed that anomalous pulse broadening, over 
a limited frequency range and possibly consis- 
tent with that observed by Lohmer et al. (2001), 
can be produced by two scattering screens hav- 
ing different scattering strengths. Such a scenario 
would produce anomalous angular broadening as 
well. I model this anomalous scattering in a piece- 
wise linear fashion, with angular broadening scal- 
ing approximately as v~ 22 below 1.4 GHz and 
above 2.7 GHz while between 1.4 and 2.7 GHz 
it would scale as v~ 1A . This model produces an 
overall frequency scaling between 0.33 and 5 GHz 
f jy-i-9 Thus, while Figure 7 shows that the an- 
gular diameter appears to scale as expected with 
frequency, a modest amount of anomalous scat- 
tering at frequencies near 1.5 GHz cannot be ex- 
cluded. 

A more stringent constraint on the anomalous 
scattering of PSR B1849+00 could be obtained 
from multi-frequency observations of its angular 
diameter. The field of view of the VLA is suffi- 
ciently large that PSR B1849+00 is detected in 
both the 0.33 and 1.5 GHz observations. Un- 
fortunately, at 0.33 GHz it is unresolved, while 
at 1.5 GHz it is sufficiently near the edge of the 
primary beam that it is difficult to obtain a reliable 
angular diameter estimate. Very long baseline in- 
tcrferometry observations are under way (and will 
be reported elsewhere) to observe B1849+005 and 
PSR B1849+00 simultaneously near 1.5 GHz in 
an effort to see if the pulsar's angular diameter 
appears anomalous. 

Nonetheless, by combining the various observa- 
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Table 5 — Continued 
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Note. — Column 5, Diameter, lists unresolved sources as .... Column 6, Structure, 
indicates whether a visual inspection of the source image suggests a contribution from 
intrinsic structure to the measured angular diameter. 



tions for B1849+005 and PSR B1849+00, it may 
be possible to constrain the properties of the re- 
gions responsible for the excess scattering. The 
scattering toward B1849+005 and PSR B1849+00 
is a factor of 5 or more in excess of what is 
predicted by the large-scale components — a thick 
disk, thin disk, and spiral arms — of the Cordes- 
Lazio model (and the model by Taylor & Cordes 
1993). In order to reproduce the observed scatter- 
ing toward PSR B1849+00, Cordes & Lazio (2003) 
introduced a "clump" or "cloud" of electrons along 
the line of sight to it. They require the clump to 
have an internal electron density n e = 10 cm~ 3 
and a fluctuation parameter F — 200. The fluctu- 
ation parameter F (see also Cordes et al. 1991) is a 
combination of the fractional variation of the elec- 
tron density internal to the clump and the outer 
scale. Larger values of F mean that a clump is 
more efficient at scattering. (For a clump, the 
fluctuation parameter and the scattering measure 
are proportional, SM ex Fn^Ar, where Ar is the 
size of the clump.) They placed the clump at a 
distance of 7.4 kpc and assumed its size to be 
Ar = 20 pc. 

The power of combining the measurements is 



twofold. First, the magnitude of angular broad- 
ening is strongly dependent upon the distribu- 
tion of scattering material along the line of sight, 
with material close to the observer contributing 
more strongly (van Langevelde et al. 1992; Lazio 
& Cordes 1998). Consider a Galactic source at a 
distance D with angular broadening 6Q a \ and an 
extragalactic source along a similar line of sight 
with angular broadening # xga i. If the scattering 
along the line of sight is dominated by a single 
screen at a (radial) distance A from the Galactic 
source, then (van Langevelde et al. 1992) 



?Gal 
'xeral 



A 
D' 



(5) 



This relationship implies that, even though seen 
through the same scattering screen, a Galactic 
source will have a smaller scattering diameter than 
an extragalactic source; this result occurs because 
the scattering screen is less effective at scattering 
the diverging wavefronts from the Galactic source 
than the planar wavefronts from the extragalactic 
source. 

The 0.33 GHz scattering diameter of B1849+005 
is 4'.' 2. Taking an upper limit on the scattering 
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diameter of PSR B1849+00 to be 3" (= 1/2 the 
beam diameter), A/D < 0.7. As PSR B1849+00 
is at a distance of roughly 8 kpc, A < 5.7 kpc 
or that the scattering screen must be more than 
2.3 kpc from the Sun. 

The second reason for combining the measure- 
ments toward B1849+005 and PSR B1849+00 
is that angular broadening and pulse broadening 
have different weightings for the distribution of 
scattering material along the line of sight. If x 
is the fractional distance between the pulsar and 
the observer (with x — being at the pulsar and 
x = 1 being at the observer), then pulse broaden- 
ing is weighted as x(l —x) while angular broaden- 
ing is weighted as x 1 . 

I have repeated the analysis of Cordes & 
Lazio (2003) but now requiring that a single 
clump describe the excess scattering toward both 
PSR B1849+00 and B1849+005. Specifically, I 
adopt the large-scale components of the Cordes- 
Lazio model and then search for a clump of excess 
scattering that can explain the observed scatter- 
ing diameter for B 1849+005, the upper limit on 
the scattering diameter of PSR B1849+00, and 
the pulse broadening of PSR B1849+00. A grid 
search over the three parameters n e , F, and x was 
used to search for the best fit in a minimum x 2 
sense. The size of the clump is assumed to be 
40 pc, so that it covers both lines of sight even 
when the clump is near the pulsar. 

An added constraint in the fitting process is 
that the dispersion measure of the pulsar (DM = 
680 pc cm~ 3 ) and the modeled distance to the 
pulsar must be self-consistent. That is, increas- 
ing the density within the clump necessarily re- 
duces the modeled distance to the pulsar, which 
in turn affects the weighting for the scattering ob- 
servables. The fitting procedure was therefore the 
following. The dispersion measure contributed by 
the clump, 5DM, was determined for each set of 
values for n e and x. The model distance for the 
pulsar was determined by requiring that the large 
scale components of the Cordes-Lazio model ac- 
count for the remainder of the pulsar's DM, i.e., 
for (DM — <5DM). This new distance was then 
used in determining the weighting factors for the 
modeled scattering observables. In practice, the 
<5DM contributed by the clump is sufficiently small 
that it does not affect the modeled pulsar distance 
(~ 9 kpc) significantly. 



No set of parameters can reproduce all three ob- 
servational constraints. In particular, the pulsar's 
pulse broadening is consistently underestimated 
(by about an order of magnitude) while the extra- 
galactic source's angular diameter is overestimated 
(by about a factor of 2-3). As might be expected, 
because the observational constraints are all due 
to scattering, the parameters n e and F are highly 
correlated. The best fit value for the distance to 
the clump is 4.7 kpc. While consistent with the 
crude estimate obtained above based on just the 
scattering diameters, this distance for the clump 
corresponds to x « 1/2. This value simply reflects 
the value of x that maximizes the modeled pulse 
broadening, which still only poorly reproduces the 
observed pulse broadening. 

If the excess scattering toward both B1849+005 
and PSR B1849+00 is attributed to the same 
clump, the angular diameters provide the most 
useful constraint; the clump must be at least 
2.3 kpc distant. Given that no single clump re- 
produces all of the observations, it is quite likely 
that different clumps affect the two lines of sight. 

Is there a structure along the line of sight 
that could contribute to anomalous scattering to- 
ward PSR B1849+00 but not B1849+005? Fig- 
ure 13 shows the Ha emission in a small region 
around these two sources taken from the South- 
ern Ha Sky Survey Atlas (SHASSA, Gaustad et 
al. 2001). Clearly apparent between the two, but 
offset slightly toward PSR B1849+00, is a clump 
of slightly stronger Ha emission. 

This Ha clump appears to lie on the outskirts 
of the H II region G33.418-0.004 (Lockman 1989). 
If this Ha clump is associated with the H II re- 
gion, the Ha clump would represent only a small 
fraction of the total Ha luminosity of the H II re- 
gion. Based on the recombination line strengths 
for the H II region, much of it must be heavily ab- 
sorbed at optical wavelengths. This line of sight 
is toward the Aquila Rift (Dame, Hartmann, & 
Thaddeus 2001), and the area to the south of the 
Ha clump appears consistent with being heavily 
obscured (Figure 13). 

The emission measure of this H II region is 
EM w 7000 pc cm- 6 (Lockman 1989). The Ha 
clump and the line of sight to PSR B1849+00 
appear to be where the 5 GHz continuum has 
dropped by approximately a factor of 2 (Altcnhoff 
et al. 1979). Thus, the EM contributed by the 
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Fig. 13. — A subimagc of Figure 12 showing the 
lines of sight to B1849+005 and PSR B1849+00. 
The gray scale shows the smoothed, continuum- 
subtracted Ha emission on a linear scale from 
to 10 Ralcighs. The sizes of the symbols are pro- 
portional to the scattering diameter (B1849+005) 
or upper limit (PSR B1849+00) at 1 GHz. For 
the scattering diameter of B1849+005, it is shown 
with the correct orientation and anisotropy 



H II region along the line of sight to the pulsar is 
perhaps 3500 pc cm~ 6 . The equivalent scattering 
measure is log SM = 0.81 [kpc m~ 20 / 3 ], assuming 
an outer scale of 1 pc (Cordes et al. 1991). This 
level of scattering is comparable to, though some- 
what below, that required for the line of sight to 
PSR B1849+00. The fitting procedure above pro- 
duces a clump with log SM = 1.4, but the scatter- 
ing parameters of the clump in the Cordes-Lazio 
model imply log SM = 3.6. 

Although it seems reasonable to attribute some 
fraction of the scattering toward PSR B1849+00 
to G33. 418— 0.004, the H II region does not appear 
to be capable of explaining it all. The line of sight 
to PSR B1849+00 would have to pierce several 
hundred similar H II regions. Higher density fila- 
ments or tendrils from H II region may be in front 
of the pulsar and responsible for the anomalous 
scattering, but the identification of such filaments 
must await higher resolution and sensitivity Ha 
observations. 

Alternately, a larger SM would be obtained 
if the outer scale in G33. 418-0. 004 were l < 
10~ 4 pc (l < 20 AU). Given that the size of 



G33. 418— 0.004 is probably comparable to 1 pc, 
it is difficult to understand why the outer scale 
would be so much smaller. Moreover, the H II re- 
gion is well offset from the location of B1849+005 
(Altcnhoff et al. 1979), so the enhanced scattering 
along the line of sight to the extragalactic source 
cannot be attributed to the H II region. As Span- 
gler et al. (1986) concluded, the enhanced scatter- 
ing toward B1849+005 (or PSR B1849+00) can- 
not be attributed to any single object or region. 

5. Conclusions 

I have observed the extragalactic source B1849+005 
at frequencies between 0.33 and 15 GHz with the 
VLA and VLBA and re-analyzed archival VLA ob- 
servations between 0.33 and 4.9 GHz. The source 
has a complex structure, as Fey et al. (1991) sug- 
gested, with components of emission both to the 
southwest and northeast of the central compo- 
nent. However, at least to 15 GHz, the central 
component is compact and highly scattered; our 
larger frequency range, nearly two orders of mag- 
nitude, extends the earlier analyses of Spangler et 
al. (1986) and Fey et al. (1991), who also found 
the source to be highly scattered. 

I have analyzed the interferometric visibilities 
to extract parameters of the spectrum of density 
fluctuations responsible for scattering and com- 
pared the scattering of B1849+005 to other nearby 
sources. The density fluctuations are anisotropic, 
with an anisotropy of about 1.3, comparable to 
that seen for other highly scattered sources in the 
Galactic disk but less than that toward Galac- 
tic center sources. While the data are sugges- 
tive that the position angle of the anisotropy 
changes with frequency, the uncertainties are 
large. Nonetheless, it seems clear that density 
fluctuations are becoming anisotropic on scales of 
order 10 15 cm (D/5kpc). An inner scale to the 
density spectrum also appears necessary to repro- 
duce the observed visibilities. The evidence for 
an inner scale l\ comes both from direct fits for a 
density spectrum including an inner scale as well 
as from assessing the values of the density spectral 
index a for a power-law-only density spectrum. 
Both of these approaches suggest an inner scale of 
roughly a few hundred kilometers, comparable to 
that found toward other sources. 

The two observations at 0.33 GHz span approx- 
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imately a decade. There is essentially no change in 
the shape of the image over this interval. A conser- 
vative upper limit on the velocity of the scattering 
material is 1800 km s _1 . This value is sufficiently 
large that it does not exclude material from the 
supernova remnant G33. 6+0.1 from contributing 
to the enhanced scattering along this line of sight. 

In the 0.33 GHz image of the field around 
B1849+005, there is a number of other sources 
that might also be heavily scattered. As is seen 
in other regions of intense scattering (e.g., the 
Cygnus region and the Galactic center), there are 
large changes in the strength of scattering on lines 
of sight separated by a degree or less. With only 
a single frequency to estimate the scattering, not 
all of the large diameters in this region can be at- 
tributed to scattering. Nonetheless, many of the 
sources probably are heavily scattered, and these 
other sources in the field could serve as a useful 
starting point for future scattering observations in 
this direction. 

Both B1849+005 and PSR B1849+00, which 
are separated by only 13', show extremely strong 
scattering, and PSR B1849+00 shows an anoma- 
lous frequency scaling for its pulse broadening be- 
tween 1.4 and 2.7 GHz. While there is no indi- 
cation of an anomalous frequency scaling for the 
angular broadening of B1849+005 above 5 GHz, 
the observations between 1.4 and 2.7 GHz (namely 
those at 1.5 and 2.3 GHz) are insufficient to place 
strong constraints on the presence of any anoma- 
lous frequency scaling. Moreover, the observations 
were insufficient for placing useful constraints on 
the angular broadening of PSR B1849+00; si- 
multaneous observations of PSR B1849+00 and 
B1849+005 at 1.6 GHz are under way and will be 
reported elsewhere. 

If the lines of sight to the pulsar and extra- 
galactic source are affected by the same "clump" 
of scattering material, it must be at least 2.3 kpc 
distant. However, a detailed attempt to account 
for all of the scattering observables toward these 
sources (angular broadening of the extragalactic 
source, pulse broadening of the pulsar, and upper 
limits on the angular broadening of the pulsar) 
does not produce a self-consistent set of parame- 
ters for a clump that can reproduce all three mea- 
sured scattering observables. 

There is a clump of Ha emission that lies in 
between the lines of sight to PSR B1849+00 and 



B1849+005, which is associated with an unob- 
scured part of the H II region G33.418-0.004. 
However, unless conditions within this H II region 
are unusual, it appears unable to produce a suf- 
ficient amount of scattering to explain the excess 
scattering along both (or even either) lines of sight. 
As Spangler et al. (1986) concluded, the scatter- 
ing toward either B1849+005 or PSR B1849+00 
cannot be attributed to a single object or region. 
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